Quantum spin ladder with dynamic phonons. Bound states of magnons and phonons. 
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We consider the two-chain Heisenberg ladder with antiferromagnetic interactions. Our approach 
is based on the description of spin excitations as triplets above a strong-coupling singlet ground 
state. It is shown that interaction with dynamic phonons drastically changes the dispersion of the 
elementary triplet excitation (magnon) and the spin structure factor. It is also shown that a new 
triplet excitation appears in the spectrum: the bound state of a magnon and a phonon. 
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Quantum spin ladders have recently become a subject 
of considerable interest, mainly because of their relevance 
to a variety of quasi one-dimensional materials |Q. It is 
well established that the two- leg S — 1/2 spin ladder 
has a gapped excitation spectrum Quite recently 

a state with the energy above the one-particle excita- 
tion level but below the two-particle continuum was ob- 
served in (VO)2P207 which can be described as a dimcr- 
ized chain 0. This result triggered theoretical work on 
two-particle bound states in this model 0. The above 
results raised the important issue of possible low-lying 
two-particle excitations in quasi one-dimensional mod- 
els with a gapped spectrum. It has been shown that in 
a two-leg ladder bound states with total spin S=0 and 
total spin S=l always exist below the continuum [^|,^). 
Similar bound states exist also in many two-dimensional 
quantum spin systems To complete the overview 

of pure spin systems (no dynamic phonons) we should 
mention that for spin chains with small explicit dimeriza- 
tion the spectrum can be even more complex consisting 
of multiple triplet and singlet states Jll],|l2| ■ 

The importance of dynamic phonons was realized a 
long time ago in relation to the spin-Pcirls transition 
p3[ . In the paper jl^| a simple model Hamiltonian for 
description of the combined spin-phonon dynamics in a 
spin chain has been suggested. Numerical study of this 
Hamiltonian has shown that many physical properties 
differ substantially from those derived in the static limit 
fusf . This important observation raises the question of 
the role of phonons in the dynamics of the spin ladder. 
In the present work, using the analytical approach de- 
veloped in |l(J, we demonstrate that dynamic phonons 
qualitatively change the spectrum of triplet excitations 
in the two-leg spin ladder: 1) an additional triplet ap- 
pears as a magnon-phonon bound state, 2) the magnon 
dispersion attains discontinuity. 3)q-dependence of the 
spin structure factor is substantially different from that 
without phonons. 

Following the idea of work pi] we consider the Hamil- 
tonian of two coupled S =1/2 chains (spin ladder) with 
local phonons coupled to the rung spin-spin interaction: 



(1) 



Here b\ and bi are the phonon creation and annihilation 
operators. We assume the optical phonon to be disper- 
sionless. The intra-chain (J) and the inter-chain (Jj_) 
interactions are assumed antiferromagnetic J, J± > 0. 

Following works ]2|,[l6| we consider the set of the rung 
singlets as the unperturbed ground state, introduce a 
triplet creation operator on the rung t* ai , a 
and hence map the Hamiltonian ([!]) to 
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where U — > oo is a fictitious interaction which enforces 
no double occupancy constraint fJrftL = 0. For shortness 
we call the excitation corresponding to v ("elementary" 
triplet) a magnon. 

We will treat the interaction with phonons by pertur- 
bation theory and therefore at the first step we have to 
diagonalize the Hamiltonian (Q) at g = 0. The technique 
to perform this diagonalization was sugg ested in Ref. §J] 
with further improvement (rainbow diagrams and attrac- 
tion in the singlet channel) developed in relation to the 
analysis of the two-dimensional J\ — J-2 model 10 1. Ap- 
plying these techniques we find that the normal Green's 
function of the magnon is of the form 



G(k,e) 



rik + is 



-e — Ofe + iS 



(3) 



Coupled Dyson's equations for the effective Bogoliubov's 
parameters Uk and Vk and for the magnon dispersion ilk 
are derived in Refs. 16 Hj. In the limit J± J the 
solution is very simple: Uk ~ 1, Vk ~ — ( J/2 J±) cos fc, 
and 



J± + Jcos k. 



(4) 



The magnon dispersion at J± = J found as a result of 
numerical solution of Dyson's equations is presented in 



1 



Fig.l (long-dashed line). For comparison we show by dots 
the results obtained from dimer series expansions 0]. 

Probing with neutrons is the experimental method 
commonly used to measure excitation spectra of mag- 
netic materials. The inelastic neutron scattering cross- 
section is directly proportional to the dynamic structure 
factor: S u (k,e) = J e l€t {S^{k, t)S%(-k, 0))dt, where we 
assume that the transverse (along the rungs) momentum 
k± = 7r, i.e. S" i = S Zt i — S' z i . Magnon contribution to 
the dynamic structure factor is S u (k, e) = S u (k)S(e— Cl k ), 
where S u {k) = (u k + v k ) 2 - At J± ^> J the structure fac- 
tor is S u (k) w 1. The plot of S u (k) at J± = J is shown 
in Fig. 2 (long-dashed line). 

To find the interaction between the magnon and the 
phonon we have to substitute transformed magnon op- 
erator t a k — Ukt a k + Vk$ a _k into the interaction term 
(g-term) of the Hamiltonian (J2J) . This gives 



Hint — ^ l-^fe.g*I,fc+ g *a. 



k,q,ct 

+ \lk, q ($ a , k+q $ a ^ k +t a ,-(k+ q )ta,k)\b q + h.c. (5) 

Tfc.g = g(uk+qUk + Vk+qVk), 
Jk.q = g(u k +qV k + U k V k +q), 

where b q is the Fourier transform of hi . 

We consider weak coupling with phonons, g 2 -C J 2 . 
Therefore the magnon self energy, given by the diagram 
in Fig. 3a, is of the form 
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Let us consider the limit J± ^> J. In this case ilk is 
given by (Q) and therefore the integration in (^|) can be 
performed analytically. This gives the correction to the 
triplet dispersion 
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This answer is satisfactory if loq > 2 J. However for lower 
phonon energy (ujq < 2 J) equation gives an infinite 
correction at k = k c , where cosfc c = (^o — J) I J- The 
reason for this is clear: the integrand in (||) has a pole and 
therefore in spite of the weak coupling the perturbation 
theory is not valid in the vicinity of k c . To improve the 
situation we have to use the exact Green's function 
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where G is defined as G(k,t) = -i{T[tk, a (t), 4,«( )]) 
jl7j . The dispersion is given by the pole of this Green's 
function. 



Note that the self energy £' in (JsJ) is different from E 
given by ([]]). The difference is very simple: to get £' 
one has to use the same eq. (^|) but with flk replaced 
by efc. This replacement is equivalent to summation of 
rainbow diagrams, Fig. 3b. Simple analysis shows that 
the vertex corrections (Fig. 3c) are negligible at small g. 
In the limit J± J eq. ([)[) has an analytical solution. 
Away from k c this solution reproduces the perturbation 
theory result (Q), however at k — k c it gives a discon- 
tinuity in dispersion. At k = k c + the dispersion is 
pushed down as e k ~ ^fe c — (g 2 /V2J) 2 ^ 3 . Immediately 
on the left from k c (k = k c — 0) the dispersion disap- 
pears because of the huge decay width. However away 
from this point (k = k c — Ak) the width drops down and 
at Ak > g A / \2J 2 Jj_ sinfc c ) the halfwidth become smaller 
than the energy. It means that the magnon reappears 
as a relatively narrow state. At J± ~ J the equation 
(^) can be solved only numerically. As an example we 
present in Fig.l (solid line) the plot of the magnon dis- 
persion at J± = J, uj = 0.5 J, g 2 / J 2 — 0.2. In spite 
of the relatively weak coupling to phonons the disper- 
sion is substantially different from that without phonons 
(dashed line). 

The magnon contribution to the the dynamic structure 
factor is sir\k)S(e- e k ), where S$T\k) = Z k {uk + v k ) 2 , 

and Zk = ^1 — ^fe-j is the quasiparticle residue. 

The plot of S^ l \k) is presented in Fig. 2 by solid 
line. It is also very much different from that without 
phonons (long-dashed line). Concluding the discussion 
of magnons we would like to note that the discontinuity 
in the spectrum is similar to the roton spectrum endpoint 
in liquid 4 He, see e.g. Ref. 

Consider now a phonon-magnon bound state. The 
binding arises due to the diagrams shown in Figs. 4a, b. 
To solve the problem in general case one has to use 
the Bethe-Salpeter equation with account of retardation. 
However for small phonon-magnon coupling (<? 2 -C J 2 ) 
which we consider in the present work, the binding is 
relatively weak and we can use simple second order 
Schrodinger perturbation theory. In this approximation 
the diagrams Fig.4a,b give the following effective inter- 
action 



M(Q,k,p) = M a (Q,k,p)+M b (Q,k,p) 



(10) 
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where Eq and Q is the energy and the momentum of the 
bound state. As soon as the interaction is known the 
bound state is defined by usual equation 

[E Q -ujo- e k ] tf(fc) = J ^M(Q, k,p)^{p). (11) 

In the limit J± ^> J the kernel is very simple M « M a « 
—g 2 /{uj + 2Jcos 2 Q/2), and solution of equation (IT 
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gives Eq = J± — J + u>o — £(,, with the wave function and 
binding energy equal to 



(8Jeg) 



1/4 



e 6 + 2Jcos 2 p/2 



(12) 



2J(w + 2Jcos 2 Q/2) 2 ' 



At J± ~ J equation (|llj) can be solved only numerically. 
We present it Fig.l (dashed line) the plot of the bound 
state dispersion at Jj_ = J, ujq = 0.5 J, g 2 /J 2 = 0.2. 
Shadowed areas show phonon-magnon continuum. It is 
clear that binding is maximum at Q = it. 

The bound state has spin S = 1 and therefore it can 
be excited in neutron scattering. The mechanism for this 
excitation is shown in Fig. 4c: first a virtual magnon is 
excited and then it emits a phonon with formation of 
the bound state. The corresponding contribution to the 
dynamic structure factor is Su\k)S(e — Ek) with Su (k) 
given by 



(Uk + Vk) 



E k -Sl k -E'(k,E k ) 



^2^k,k- P ip(p) 
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At Jx 3> J integration in this formula with account of 
eqs. (Oh gives 



Si b \k) 



(w -2 J cos 2 k/2f 
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It is assumed that the denominator in this formula is not 
too small. At small denominator there is strong mixing 
between the magnon and the bound state and therefore a 
more accurate consideration is necessary. The structure 
factor S ( u b \k) for J ± = J, u Q = 0.5 J and g 2 = 0.2 J 2 
obtained by numerical integration in ([l3]) is shown in 
Fig. 2 by the dashed line. At k = it the bound state 
excitation probability is 10 times smaller than that for 
the magnon. However at k ~ 0.87T the probabilities are 
comparable. 

Above we have considered the magnon and the bound 
state of the magnon with a phonon. In addition to these 
states there are bound states of two magnons considered 
in Ref. ||. The dynamic phonon does not qualitatively 
influence these bound states. However additional phonon 
diagrams shown in Fig. 5 push down the energies of these 
states. To illustrate this for J± — J, to — 0.5J and 
g 2 = 0.2J 2 we show in Fig.l (dashed-dotted line) the 
dispersion of the two-magnon bound state with S = 1. 
Because of the decay to the phonon and two magnons 
this state attains finite width T ~ 0.15 J This 
is why the state disappears at k ~ 0.6-7T where split- 
ting from the corresponding continuum is getting smaller 
than halfwidth. The structure factor for this state is 
S g (k,e) = S g (k)Sr(e — Ek), where Sr is a function with 



finite width. The index g indicates that the state can be 
excited at k± = 0, i.e. by Si + S-. The formula for S g (k) 
is presented in ||. Integration in this formula gives the 
plot shown in Fig. 2 by the dotted-dashed line. At k = it 
the excitation probability is 4 times smaller than that for 
the magnon. However at k ~ 0.8vr they become equal. 

In conclusion, we have studied the two-leg antiferro- 
magnetic Heisenberg ladder with spins coupled to dy- 
namic phonons. In the analysis we use the bond operator 
representation together with the diagram technique and 
Brueckner approximation. It has been demonstrated that 
the phonons qualitatively change spectrum of triplet ex- 
citations: 1) An additional triplet appears as a magnon- 
phonon bound state. Dispersion of this triplet for some 
particular parameters is shown in Fig.l (dashed line); 2) 
the magnon (elementary triplet) dispersion attains a dis- 
continuity, Fig.l, solid line; 3)k-dependence of the spin 
structure factor (Fig. 2, solid line) is substantially differ- 
ent from that without phonons (Fig. 2, long-dashed line). 

I would like to thank M. Kuchiev for stimulating dis- 
cussions. 
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FIG. 1. Triplet excitations on the isotropic (J± — J) 
ladder. The long-dashed line is the magnon dispersion ob- 
tained by the Brueckner method without coupling to phonons. 
Dots represent the same dispersion obtained in Ref. [5] us- 
ing series expansion. Solid line gives the magnon dispersion 
with account of coupling to phonons. The phonon energy 
is wo = 0.5J and the coupling constant g 2 /J 2 = 0.2. The 
dashed line corresponds to the bound state of the magnon 
and the phonon. The dashed-dotted line corresponds to the 
bound state of the two magnons with total spin S=l. Shad- 
owed areas gives magnon-phonon continuum. 

FIG. 2. Spin structure factors (excitation probabilities) for 
isotropic ladder, J± — J. The long-dashed line corresponds 
to the magnon without coupling to phonons. Solid line de- 
scribes the same magnon coupled to phonons (o»q = 0.5J, 
g 2 / J 2 = 0.2). The dashed line describes bound state of the 
magnon and the phonon. The dashed-dotted line corresponds 
to the bound state of the two magnons. 

FIG. 3. Magnon self energy due to the phonons. Dashed 
line corresponds to the phonon Green's function, (a) Single 
loop, (b) rainbow diagrams, (c) vertex correction. 

FIG. 4. Diagrams (a) and (b) represent the phonon - 
magnon scattering amplitudes responsible for the binding. 
Solid line describes the magnon and dashed line describes the 
phonon. The diagram (c) shows mechanism for the bound 
state contribution to the spin structure factor. 

FIG. 5. Phonon corrections for the two magnon bound 
state. 
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